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Abstract 

We analyze a simple dynamical model of glasses, based on the idea that each 
particle is trapped in a local potential well, which itself evolves due to hopping 
of neighbouring particles. The glass transition is signalled by the fact that 
the equilibrium distribution ceases to be normalisable, and dynamics becomes 
non-stationary. We generically find stretching of the correlation function at 
low temperatures and a Vogel-Fulcher like behaviour of the terminal time. 
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Glasses have a number of fascinatingly universal properties which are still not satisfacto- 
rily accounted for theoretically [I], 0. A common experimental feature is the 'shouldering' 
of the relaxation laws. More precisely, the relaxation of - say - the density fluctuations 
evolves from a simple Debye exponential at high temperatures (liquid) to a two-step process 
at lower temperature, where the correlation function decays fast to a plateau value, from 
which it subsequently decays on a much longer time scale. These two regimes are called, 
respectively, the (3 and a relaxations; the a decay is often described in terms of a 'stretched 
exponential' with a characteristic time scale r diverging faster than exponentially with the 
temperature, and controlling the transport properties such as the viscosity. The most pop- 
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ular description of this divergence is the Vogel-Fulcher law: r ~ r e T ~ T o, where r is a 
microscopic time scale |J. One has to note that despite its tremendous phenomenological 
success, this law predicts such an abrupt divergence when T is lowered that it cannot be 
tested near Tq. Correspondingly, other functional forms, such as r ~ Toe^-* 2 , give reasonable 
fits of the data 0,0. Furthermore, the Vogel-Fulcher law has only been justified on rather 
heuristic grounds ||. 

Up to now, the most comprehensive theory of dynamical processes in glasses is the so- 
called mode-coupling theory, developped by Gotze and others ||. It is based on a family 
of schematic equations coupling the density fluctuations in a non-linear and retarded way. 
Generically, these equations have a singularity which is associated to an 'ideal glass' transi- 
tion temperature T c , below which the correlation function does not decay to zero ('broken 
ergodicity'). This theory describes satisfactorily the overall shape of the relaxation function, 
at least for T > T c - in particular the existence of the two regimes and a mentionned 
above, and a power-law divergence of the 'terminal' time scale r as (T — T^) -7 . How- 
ever, comparison with experiments |7| shows that the transition temperature T c , if it exists, 
is much higher than the Vogel-Fulcher temperature, leaving a whole temperature interval 
[To,T c ] where mode-coupling predicts a 'fluctuation arrest' (i.e. the a regime disappears) 
while the experimental relaxation time is still finite (and behaves a la Vogel-Fulcher). A 
way out of this contradiction is to argue that the mode coupling theory leaves out 'acti- 
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vated processes' which are responsible for the long time relaxation, and act to blur out the 
power-law divergence of r near T c . Although not unconceivable, this possibility requires the 
introduction of at least one extra free parameter to fit the data, which would be zero in the 
ideal transition scenario and is not found to be small in the experiments. In other words, one 
major aspect of the mode-coupling theory is the existence of a singular temperature which 
however does not manifest itself very directly experimentally - in particular, the terminal 
time r does not reveal any accident around T c [0. 

Another - rather more subtle - difficulty associated with the mode-coupling theory is that 
the dynamical equations are formally identical ||, |§ to those describing exactly some mean 
field models of spin-glasses ||, [[UJ, where the presence of quenched disorder is assumed from 
the start. In glasses, however, this quenched disorder must be in some sense 'self-induced'. 



Although some progress has recently been made to substantiate such a scenario [pj] , [ |12|j , it 
is not yet clear whether the glassiness found in mode-coupling theories is or not an artefact 
of the very approximation. 

In this paper, we propose and solve a simple model of glasses. Although still rather ab- 
stract, we believe that it captures at least part of the physics involved in the glass transition. 
The shape of the correlation function evolves, as the temperature is decreased, precisely as 
in experimental glasses - in particular, the terminal time diverges according to the Vogel- 
Fulcher law. The glass transition is signalled by the fact that the equilibrium distribution 



ceases to be normalisable; correlatively, as argued in |L3|], [p~4|| , aging effects are present in 
the glass phase. 

The progressive freezing of a liquid can be thought as follows: each particle is in a 
'cage', i.e. a potential well of depth e created by its neighbours, from which it can escape 
through thermal activation. However, since a priori all particles can move, the (random) 
potential well trapping any one of them is in fact time dependent, further enhancing the 
probability of moving. In order to understand the glass transition, one must describe how, 
in a self-consistent way, all motion ceases. We thus introduce a density of local potential 
depth p(e), describing the fact that the efficiency of the 'traps' depend on the environment 



|15|| . Now, the basic object on which we shall focus is the probability P(e,t) that a given 
particle is in a trap of depth e at time t. This probability evolves because a given particle, 
with rate r exp — ^, leaves its trap and chooses a new one with weight p(e). Doing so, all 
the neighbouring 'traps' are affected by the motion which has taken place. In a mean field 
description, the resulting evolution of P(e, t) is described by the following equation: 



-r„exp-(l) P(e, t) + T(t)p(e) + !•(*)£>!■ 



p(e) ^M_P (M) .^) 



de de 
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dt " r T v " wrw w & 

where T(t) = ro(exp— (^)) is the average hopping rate ((...) means an average over P{e,t) 
itself). The two first terms describe the direct effect of leaving a trap, while the third 
one expresses the fact that every 'hop' induces a small change in all the neighbouring e's. 
Assuming that the transition rate is proportional to the density of final states, the balance 
equation reads: T(t) J de' T(e — e'){P(e', t)p(e) — P(e, t)p(e')}. The fact that the change is 
small, justified in a mean-field limit where the number of neighbours is large, allows one to 
write this term in a diffusion like fashion, with an effective diffusion constant D proportional 
to the width of T. With no incidence on the following results, we shall restrict e to be positive 
- in line with our trap picture. Eq. (1) is then supplemented by the boundary condition: 

(,(e)^-P(M)^)=0 for £ = (2) 

which means that e = is a 'reflecting' point. In fact, Eqs. (1) and (2) can be taken as 
a definition of our dynamical model for glasses, which must be supplemented by an initial 
condition P(e, t = 0). 

Immediate properties of Eqs. (1,2) are that: 

• J-m deP(e, t) is a conserved quantity, as it should. 

• When T — > oo, exp — (^) = 1 and the equilibrium distribution is simply given by 
P eq (e) = p(e), as expected since there is no Boltzmann factor biaising the a priori weights. 

Let us now study the equilibrium distribution at finite T. Setting 9Pc Q t ^ = 0, one obtains 
an inhomogoneous Schrodinger equation for P eq {e): 
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_ D ^(£) + V( £ )P« q (e)^ v( ^ r„ex P -(g.r^ (3) 

with the boundary condition Eq. (2). As long as this equation admits a 'bound state', 
a normalisable -P e q(e) exists and we shall call the resulting state 'liquid'. However, as the 
temperature decreases, the effective potential V tends to push P eq (e) towards larger e, and 
depending on the shape of p(e), an 'extended', non normalisable state may appear - corre- 
sponding to a glass phase. It is easy to show that if p(e) decreases slower than exponentially 
for large e, the bound state ceases to exist as soon as T < oo, while if p(e) decays faster than 
exponentially, the bound state remains down to T = (we shall come back to this case be- 
low). We shall thus focus on the case where p(e) is a simple exponential: p(e) = ^ exp —jr, 
where T turns out to be the glass transition temperature. The reason for this is quite 
simple: it is the temperature at which the Boltzmann weighting factor exactly compensates 
the fact that deep potential wells are a priori extremely rare. Such a scenario is reminiscent 
of Derrida's random energy model [l6j , where the transition temperature is also defined by 
balancing the locally exponential density of states with the Boltzmann factor. 

Solving Eq. (3) for -P e q( e ) m this case leads to a well defined equilibrium state for T > T , 
which reads fl7fl: 



^eq(e) 



J\f 



/ •) !'.' £„(■>•„) + K v {x){X u {x Q ) -X v {x)) + I v {x)K, v (x) 



K v _i(xq) 



(4) 



where v 



2T 



n 3/2 



X 



Xoexp-^r, and xq = uT %^ :y ^. I u and K v are the Bessel functions 



t-t > - — -u —i- vTo , " V5 Vf 

roo du j/- („,\ t ( ™\ — r x du 



of order u, and K v (x) = f K v (u), Z v {x) = f* fl„(u). M and V are fixed by the 
normalisation of -P e q(e) and the boundary condition (2), which lead to the following equation: 

D I v -i{x Q ) o _ f°° du 



^-[Mx )] 2 + 2 r —IMK v {u) 

l{Xo) Jx U 



(5) 



One can check, using the properties of Bessel functions, that this last equation is identically 
satisfied when T — > oo, where v = 2 and r = Tq. For T — > To, on the other hand, we find 
that the average hopping rate T vanishes linearly, as r T ~ T ° , which is numerically found to 
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be a very good approximation for all temperatures. More interestingly, however, one finds 

T ~ T ° r-n £(T ~ T( 
TT ^t> TT Q 



that for T — > T , -P eq (e) decays as P eq (e) ~ ^ff 3 ' ex P ~ ^tt 0' which means that the 



characteristic energy scale is e* = ■ 

In order to make contact with experimental observables one must further define a (two- 
time) correlation function. The simplest one to consider, corresponding to a high dimensional 
space, is such that only particles which have not moved at all contribute to the correlation, 



i.e. 



poo f 

C(t,t') = J o deP(e,t')exp-[T exp-( f )(t-t')] (6) 

but other choices, corresponding e.g. to particles hopping on a d-dimensional lattice, are 
possible | [2 Of . Eq. (6) assumes in particular that the 'width' of the potential wells is zero. In 
order to be more realistic and take into account the fast vibration of the particles in their 
'cages', a simple modification is to multiply C(t,t') defined in Eq. (6) by a 'Debye- Waller' 
factor Cp(t, f) = exp — 9 r ( ^~ t - , where q is the probing wave vector. r(t) describes a diffusive 
motion in an harmonic potential well: r 2 (t) = £q[1 — exp(— ^)]; £o can be thought as the 
'size' of the cage, and of the order of the high temperature (liquid) diffusion constant |21 . 

For T > Tq, i.e. when P e q{e) exists, the correlation function only depends on the differ- 
ence t — t'. One finds that C(t) (defined by Eq. (6)) behaves as 

r i - rt for t < r- 1 

C(t) = ^ (7) 

[(V t) T « for t>r(T) 

where r(T) = Tq 1 exp( T ^° To ) is the Vogel-Fulcher time, which very naturally appears within 
the present model (although no diverging length scale is involved). From Eq. (7), and 
Fig. 1, one sees that C(t)Cp(t) has precisely the shape observed in most experimental 
situations, provided one takes into account the Debye- Waller factor Cp defined above. Note 
the presence of two characteristic time scales, a short (microscopic) one To, corresponding 
to the cage vibrations (ft peak), and a long one r(T), corresponding to the a peak; these 
two time scales separate extremely fast as the temperature is reduced. 

When T < Tq, on the other hand, no normalisable -Peq(e) can be found, which corresponds 
to the weak ergodicity breaking situation described in [[Tj|, [|14], flO l- In this situation, 
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P{e,t) never reaches a stationary limit, but continuously drifts towards larger and larger 
energies. Time translational invariance is spontaneously broken as C(t,t') never becomes a 



function of t — t' alone, a situation now referred to as 'aging'. In fact, one can show (18 
that C(t w + t,t w ) = C{ — ), with 1 — C{u) oc u T o for u — > and C oc u T o for u — > oo, 
precisely as in the 'trap' model studied in ||13|| . This suggests that the difference between the 
'quenched' model considered in [|13| (corresponding to D = in Eq. (1)) and the 'annealed' 
model considered here is, to some extent, irrelevant. 

All the above results are still expected to hold if the density of states p(e) is approx- 
imatively exponential below a certain cut-off e c , provided that T — To > If p(e) 

decays faster than exponentially, for example as p(e) oc exp — {j-) 2 , then strictly speaking 
the glass temperature is pushed down to zero. This choice for p(e) is interesting since it 
really corresponds to a mean-field limit where the local trap strength is obtained as a sum 
of contributions from the (large) number of neighbours. The equilibrium distribution -P e q(e) 
is then, for large e, given by exp ^ — (j-) 2 , and the corresponding C(t) exhibits a con- 
siderable amount of stretching at low enough temperatures. C(t) is indeed very well fitted 
by a stretched exponential at intermediate times ||22|| . The long time fall off of C(t) is in 
this case given by C(t) oc {^p-)* 1 with fi = (^) 2 logr t. The terminal time r(T) diverges 
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as exp|%, i.e. much faster than an activated law, and, as mentioned in the introduction, 
also compatible with the experimental data [Q, ||. The shape of C(t)Cp(t) is plotted in 
Fig. 2 for different — . It is, again, very similar to the experimental data. In particular, the 
a relaxation for different temperatures can be approximatively |23[ rescaled onto a unique 



master curve when plotted as a function of t^j, as observed in experiments and numerical 
simulations [24[, and predicted by the mode-coupling theory 0. 



In conclusion, we have proposed a very idealized model for the glass transition, which we 
have solved in the high temperature phase. The results are found to reproduce two striking 
observations: the shouldering and stretching of the correlation and the Vogel-Fulcher like 
divergence of the terminal time scale. The former feature is usually accounted for by the 
mode-coupling theory on the basis of a 'phantom' singularity which is removed by exogeneous 



processes, in turn responsible for the fmiteness of the terminal time at low temperatures. 
Our model suggests that the hypothesis of a intermediate temperature transition might 
not be necessary, although a more detailed comparison with experimental data is obviously 
desirable. 
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Cugliandolo, D. Dean, J. Kurchan, D. Lairez, J. M. Luck, E. Vincent and especially M. 
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Figure Captions. 

Fig 1: Exponential density of states. Plot of C(t)Cp(t) versus log 10 (r t) for q£ = 0.5, 
r = 5To \ and § = 2., 1.1, 1.03. 

Fig 2: Gaussian density of states. Plot of C(t)Cp(t) versus log 10 (r t) for q£ = 0.5, t = 
5Fq 1 , and j- = 0.5, 0.25, 0.15. Note that the plateau observed for the lowest temperature 
eventually decays to zero. 



8 



REFERENCES 



[1] For a review, see the interesting series of papers in Science, 267 (1995) 1924. 

[2] 'Liquids, freezing and glass transition', Les Houches 1989, JP Hansen, D. Levesque, J. 
Zinn- Justin Editors, North Holland, in particular ref. H 

[3] H. Vogel, Z. Phys. 22 (1921) 645 ; G. S. Fulcher, J. Am. Ceram. Soc, 6 (1925) 339 

[4] J.D. Ferry, L.D. Grandine, E.R. Fitzgerald, J. Appl. Phys. 24 (1953) 911, H. Bassler, 
Phys. Rev. Lett 58 (1987) 767 and refs. therein. 

[5] W. Gotze, in ref. [2], p. 403 et sq. 

[6] see T. R. Kirkpatrick, D. Thirumalai, P. G. Wolynes in ref. ||; T. Vilgis, J. Phys. Cond. 
Matter, 2 (1990) 3667, G. Parisi, Slow dynamics in glasses, preprint cond-mat 941115 
and 9412034. 

[7] for a enlightening introduction to the experimental controversy, see the series of Com- 
ments in Phys. Rev. E: X.C. Zeng, D. Kivelson, G. Tarjus, Phys. Rev. E 50 (1994) 1711, 
P. K. Dixon, N. Menon, S. R. Nagel, Phys. Rev. E 50 (1994) 1717, H. Z. Cummins, G. 
Li, Phys. Rev. E 50 (1994) 1720, and references therein, in particular H. Z. Cummins, 
W.M. Du, M. Fuchs, W. Gotze, S. Hildebrand, A. Latz, G. Li, N.J. Tao, Phys. Rev. E 
47 (1993) 4223 

[8] T. R. Kirkpatrick, D. Thirumalai, P. G. Wolynes, Phys. Rev. A 40 (1989) 1045 and 
references therein. See also Th. Nieuwenhuizen, preprint (1995). 

[9] J. P. Bouchaud, L. Cugliandolo, J. Kurchan, M. Mezard, in preparation. 

[10] L. Cugliandolo, J. Kurchan, Phys. Rev. Lett. 71 (1993) 173, L. Cugliandolo, P. Le 
Doussal, Large time off- equilibrium dynamics of a particle diffusing in a random poten- 
tial, preprint |cond-mat 95051 12| . S. Franz, M. Mezard, Europhys. Lett. 26 (1994) 209, 



Physica A209 (1994) 1 



9 



[11] T. R. Kirkpatrick, D. Thirumalai, J. Phys. A 22 (1989) L149. 

[12] J. P. Bouchaud, M. Mezard J. Phys. I France 4 (1994) 1109, E. Marinari, G. Parisi, and 
F. Ritort, J. Phys. A 27 (1994) 7615 ; J. Phys. A 27 (1994) 7647 , L. F. Cugliandolo, 
J. Kurchan, G. Parisi, F.Ritort, Phys. Rev. Lett. 74 (1995) 1012 , S. Franz, J. Hertz, 
Phys. Rev. Lett. 74 (1995) 2114 . See also: F. Ritort, Glassiness in a Model without 
energy barriers, preprint cond-mat 9504081 

[13] J. P. Bouchaud, J. Physique I 2 (1992) 1705, J. P. Bouchaud and D. S. Dean, J. Phys. 
I (France) 5 (1995) 265 

[14] F. Bardou, J. P. Bouchaud, O. Emile, C. Cohen- Tannouji, A. Aspect, Phys. Rev. Lett. 
72 (1994) 203 

[15] Note that this one-particle density of state can a priori depend on the density of the 
liquid and the temperature itself. 

[16] B. Derrida, Phys. Rev. B 24, (1981) 2613 

[17] The problem of the convergence towards this equilibrium state is interesting and will 



be discussed in [ITS 



[18] J. P. Bouchaud, A. Comtet, C. Monthus, in preparation. 

[19] This simple exponential form is actually a good numerical approximation for all values 
of e. 

[20] see e.g. J.P Bouchaud, A. Georges, Phys. Rep. 195 (1990) p. 142-143. 

[21] This model emerged from discussions with M. Lairez and M. Adam: see D. Lairez, E. 
Raspaud, M. Adam, J.P. Carton, J.P. Bouchaud, to appear in 'Die Makromol. Chem. 
Symp.' (1995). 

[22] On this point, see, e.g. B. Castaing, J. Souletie, J. Phys. I (France) 1 (1991) 403. 



10 



[23] Note that fi(t = r(T)) is indeed independent of temperature. 
[24] W. Kob, H.C. Andersen, Phys. Rev. Lett. 73 (1994) 1376 



11 



